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Abstract Let A be an {m x n) integral matrix, and let P = {x : Ax < b} 
be an n-dimensional polytope. The width of P is defined as w{P) = min{x € 
Z"\{0} : maxx^px^u — minx^px^v}. Let A{A) and 5(A) denote the greatest 
and the smallest absolute values of a determinant among all r{A) x r{A) sub¬ 
matrices of A, where r(A) is the rank of a matrix A. 

We prove that if every r(A) x r(A) sub-matrix of A has a determinant 
equal to ±Z\(A) or 0 and w{P) > (A(A) — l)(n -I- 1), then P contains n 
affine independent integer points. Also we have similar results for the case of 
k-modular matrices. The matrix A is called totally k-modular if every square 
sub-matrix of A has a determinant in the set {0, ±fc’’ : r G N}. 

When P is a simplex and w{P) > 5(A) — 1, we describe a polynomial time 
algorithm for finding an integer point in P. 

Finally we show that if A is almost unimodular, then integer program 
maxjc^x : x € P Ci Z"} can be solved in polynomial time. The matrix A 
is called almost unimodular if Z\(A) < 2 and any (r(A) — 1) x (r(A) — 1) 
sub-matrix has a determinant from the set {0, ±1}. 
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1 Introduction 

Let A be an m X n integral matrix. Its ij-th element is denoted by Aij, 
is the i-th row of A, and is the j-th column of A. For a vector 6 G Z”, by 
P{A, b) (or by P) we denote the polyhedron {x € K" : Ax < b}. The set of all 
vertices of a polyhedron P is denoted by vert{P). 

Let r{A) be the rank of a matrix A. Let A{A) and S{A) denote the greatest 
and the smallest absolute values of the determinant among all r{A) x r{A) sub¬ 
matrices of A. Let Aicm{A) and Agcd{A) be the least common multiple and 
greatest common divisor (resp.) of absolute values of a determinant among all 
r{A) X r{A) sub-matrices of A. 

Now we refer to the notion of the k-modular matrices that has been in¬ 
troduced in m- The matrix A is called totally fc-modular if every square 
sub-matrix of A has a determinant in the set {0, ifc’’ : r G N}. The matrix A 
is called k-modular if every r{A) x r{A) sub-matrix of A has a determinant in 
the set {0, ifc’’ : r G N}. 

Also we refer to the notion of almost unimodular matrices that was in¬ 
troduced in [5] for square case. The matrix A is called almost unimodular if 
A{A) < 2 and any (r(A) — 1) x {r{A) — 1) sub-matrix has a determinant from 
the set {0, ±1} 

For a matrix B G cone{B) = {x : x = Bt, t G K", ti > 0} is a 

cone spanned by columns of B and conv{B) = {x ■. x = Bt, t G M", U > 
0, X]r=i ~ convex hull spanned by columns of B. 

For a vertex v of P{A,b), N{v) = {x G R." : Aj t^x < b}, where J = {k : 
Ak*v = b}. 

The following theorem was proved in |35j 

Theorem 1 If every n x n determinant of A belongs to {—2, —1, 0,1, 2} and 
P{A, b) is full-dimensional, then 

1. P(A,&)nZ” ^0. 

2. One can cheek the emptiness of the set P{A, b)r\JP in polynomial time. 

3. For every row a of A, the problem max{a^x : Ax < b, x € Z"} can be 
solved in polynomial time. 

4 . For every v G vert(conv{P{A,b) DIP)) there exists u G vert{P), such 
that V lies on some edge which contain u (v lies on an edge of N(u)). 

5. If each nxn sub-determinant of A is not equal to zero, then the problem 
max{c^X : Ax < b, x € Z"} can be solved in polynomial time. 

The remarkable result was obtained by V.E. Alekseev and D.V. Zakharova in 
[U for {0, l}-matrices. 

Theorem 2 Let A G {0,1}"^^", b G {0,1}™, c G {0,1}". Let all rows of A 
have at most 2 ones. Then, the problem max{c^x : Ax < b, x G Z"} can be 
solved in polynomial time when Z\(( ^ }) is fixed. 

Width of P is defined as w{P) = min{x G Z” \ {0} : maXx^px^u — 
minx^px^ v}. 
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Now we refer to the classical flatness theorem due to Khinchine m- Let 
P be a convex body. Khinchine shows that if P fl Z" = 0, then w{P) < 
f{n), where f{n) is a value that depends only on a dimension. There are 
many estimates on /(n) in the works [HISlFTl fTTlI^IST] . There is a conjecture 
claiming that f{n) = 0{n) [71IT7]. The best known upper bound on f[n) is 
log‘^{n)) due to Rudelson [5T] . where c is some constant that does not 
depend on n. 

The paper m contains an estimate of the width for a special class of 
polytopes. 

Theorem 3 Let A € b € Z™, P{A,b) is a polytope and every r{A) x 

r{A) sub-determinant of matrix A is equal to ■iiA(A) or 0. If w{P{A,b)) > 
{A{A) — 1) (n + 1), then \P{A,b) fl Z"| > n + 1. Moreover we can find an 
integer point in P{A, b) H Z" using a polynomial time algorithm. 

We give here another proof of this result. 

An interesting problem is estimating of f(n) for empty lattice simplices [H 
[T3l[T8l[32] . A simplex S is called empty lattice if vert{S) C IP and S' n Z” \ 
vert{S) = 0. Best known estimate of f{n) for the empty lattice simplices is 
0 {nlog{n)) due to [3]. 

In this paper we will prove that the width of a simplex (not necessary with 
integer vertices) without lattice points is at most 5{A) — 1, where A is the 
restriction matrix of the simplex. Moreover, if its width is at least (5(A) —1, then 
we can find an integer point in the simplex by a polynomial-time algorithm 
presented in this paper. 

The authors consider this paper as a part of general problem for finding 
out critical values of parameters, when a given problem changes complexity. 
For example, the integer programming problem is polynomial time solvable on 
polyhedrons with integer vertices, due to m- On the other hand, it is NP- 
complete in the class of polyhedrons with denominators of extreme points equal 
1 or 2, due m- The famous fc-satisfiability problem is polynomial for k < 2 
but is NP-complete for all fc > 2. In the papers some graph parameters 

(the density and packing number) were considered and it was described how 
its growth in terms of the number of vertices affects on the complexity of 
the independent set problem. A theory, when an NP-complete graph problem 
becames easier, is investingated appling to the family of hereditary classes 
in the papers Our main interest is to determine a 

dependence of the integer programming problem complexity on spectrum of 
sub-determinants of the restriction matrix. 


2 The polytopes with bounded determinants 


The main result of this section is 
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Theorem 4 Let A G b G Z"* and P = P{A, b) is a polytope. Ifw(P) > 


{Ai,^{A) - !)■ 


Z\(A) 


■{n + 1), then \P{A, b) fl Z"| > n + 1. We can find an 


AgcdiA) 

integer point in P{A, b) H Z" using a polynomial time algorithm. 


This theorem has two trivial corollaries: 


Corollary 1 Let A G Z™^" be a fc-modular matrix, b G Z™ and P = P{A, b) 
is a polytope. 

If w{P) > {A{A) — 1) + 1), then \P{A, 6) n Z"| > n + 1. We can 

d{A) 

find an integer point in P{A, b) n Z" using a polynomial time algorithm. 

Corollary 2 Let A G Z™^", b G Z™, P = P{A,b) is a polytope and each 
r(A) X r(A) sub-determinant of A is equal to ±A(A) or 0. 

If w{P) > (A{A) — l)(n + 1), then \P{A, b) n 1A\ > n + 1. We can find an 
integer point in P{A, b) n Z” using a polynomial time algorithm. 

We need the following lemmas to prove the theorem |4l 

Lemma 1 Let B € , M = cone{B)r\'IP\{0}, and x* = minx^MC^x, 

then X* G conv{0B). 

Proof Suppose, that c^B^,i = miUiC^B^,i and x* = Bt ^ conv(0,bi, ...,bs), 

then ti H- \-ts > 1. Therefore x* = Bt > (ti H- B.,i, 

this is a contradiction. 

Lemma 2 Let z' = b'"^x' = min{b'^x : x G M}, z" = b"^x" = min{b"^x : 
X G M}, then \z' — z"\ < \\b' — 6"||oo||a:”||i- 

Proof Let z' > z", then z' - z" < b'^x” - b"^x” < ||6' - &"||oo||a:"||i. 
Lemma 3 Let P' = P{A, b') and P" = P{A, b") are nonempty polytopes, 
then \w{P') - w{P'')\ < + 1)11^' “ ^"lU • 

Proof w{P') = min{maxp'C^x — minpicAx : c G Z” \ {0}}. According to the 
duality theorem of the linear programming 

maxpic^x = min{b'^y : A^y = c, j/ > 0} 

minpic^X = —min{b'^z : —A^z = c, z > 0}. 

Therefore, 

w{P') = min{b'^y + b'^z : A~^y + A^z = 0, y > 0, z > 0, A^y G Z" \ {0}}. 

Let’s consider the cone C = {(^^ G : A^y + A^z = 0, (^) > 0}. From the 
equation A^y + A^z = 0 and the Kramer’s rule it follows that C is generated 
by the vectors with components in the set of all r(A) x r(A) sub-determinants 
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of A. Therefore, the maximal absolute value of coordinates of these vectors is at 
most A{A). This value can be decreased by a dividing each of components by 
Agcd{A). According to Carateodory’s theorem, every vector has at most n + 1 
non-zero coordinates. So, by the lemma [TJ there are y',z' such that w{P') = 

b'^y' + b'^z' and -|-Z-) < ——Similarly we can show that 

^gcd[A) 

w{P") = b"^y" + b"^z" and YZi iv” + 4') < + 1) for some y", z". 

‘^gcd(A) 

Suppose, for clearness, that w{P') > w{P"). We have \w{P') — w{P”)\ = 
{b'^y' + b'^z') - b"^z") < {b'^y" + b'^z") - ib"^y” + b"^z") = 

(b'-h"yy"+{b'-b"yz" < ii6'-6"iiooii(y",z")lli < 

Now, we are ready to prove the main result (theorem 0]) of this section. 


Proof Let b' e Z"*, such that b'^ = bi — {bi mod Aicm{A)) and P' = P{A, b'). 
If w{P) > (A[cm{A) — 1)-;—-|- 1), then by the lemma |31 w{P') > 0. 

^gcd[A) 

Thus, P' is full-dimensional and each component of b' is divided by Aicm{A). 
So, it is easy to see, that P' is a full dimensional polytope and all components 
of any vertex of P' are integer. Now the theorem follows from the fact that 
P' C P. 

We can use any polynomial algorithm of linear programming (Khachiyan’s 
algorithm [ISKS]) to find some vertex of P' as an integer point of P. 


3 The simplices with bounded determinants 

A part of this section describes results of the R. E. Gomory [SlITUligilTl] . We 
will repeat some of the Gomory’s arguments, slightly modifying them for our 
purposes. 

Let A € B e Z”><*, b e Z", and \det{A)\ = A > 0. Gonsider the 

system 


\Ax + By = b 

jccGZ’^, yeZ^ ■ ^ ^ 

Let D be the Smith normal form [53] of the matrix A, then A = P~^DQ~^, 
where P~^, Q~^ are integer unimodular matrices. So system © becomes 


DQ + PBy = Pb 
x G Z", 2/ G Z^ 
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After the unimodular map Q ^ x and removing x’s variables the sys¬ 
tem becomes 


\ PBy = Pb (mod D) 

There is a bijection between variables x and y giving by the formula x = 
A-\b-By). 

Let M{A, B, b) be a polyhedron induced by the system ([2]). 

Definition 1 (Gomory) Let y G M{A, B, b). Hence PBy = Pb(modD). 
We say that y is an irreducible point of M(A, _B, b) if for any u ^ v such that 
u < y and v < y we have PBu ^ PBv (mod D). 

Theorem 5 (Gomory) Columns of the matrix PB induce an additive group, 
the group operation is an addition by a modulo D. The power of this group is 
at most A. If M(A, B, b) ^ (h and y is an irreducible point of M(A, B, b), 
then nfci(l + Vk) < 

Proof It is easy to see that columns of PB induce an additive group. Let 
g be element of this group, so g = PBt mod D for some t G Hence 
0 < gk ^ Dkk — Ij where 1 < fc < n. So a total number of group elements is 
at most nLi^fefc = ^- 

Let y be an irreducible point of M(A, B, b), and t < y. From the definition 
of the point y it follows that all group elements g = PBt mod A are different 
for different t, so different combinations of t induce different elements of a 
group. Since the total number of t combinations is nLi(l + Vk) and the 
number of distinct group elements is at most A then we have 0^=1 (1 + Vk) ^ 
Z\. 

Theorem 6 (Gomory) Let y G vert(M(A, B, b)), then y is an irreducible 
point of M (A, B, b). 

Gorollary 3 (Gomory) Let y G vert(M(A, B, b)), then 0^=1 (1 +J/fc) < 

Lemma 4 Let I be an identity matrix. There is polynomial time algorithm to 
find the point y G M(A, I, b) with properties< A — 1 andYly^^yoPk < 
- 1 . 

Proof The system for M (A, I, b) has very simple structure: 


Py = Pb (modD) 

yGl% 
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So y = b + P '^Dt, for any t G If. We can assume that the matrix P 
l^sxs following form; 


-1 



P2\ 1 

0 ... 

0 

p-l = 

PSl PS 2 

1 ... 

0 


\Psl Ps2 Ps3 ■ ■ ■ 

y 


Hence 2/1 = and j/fc = bk + {pkiDiiti+pk 2 D 22 t 2 -\ - \-Dkktk) = 

bk+Dkktk + z2j=i PkjDj jtj for 1 < fc < s. Then we can put yi = bi mod Di i 
and yk = {bk + PkjDj jtj) mod ^ for 1 < fc < s. 

Hence we have vector y G M{A, /, b) with the property OLi Vk OLi 
1). We know that -1 Dk k = A and a maximum of the sum ^kk, that 

equals to Z\ + s — 1 is reached when Dk fe = 1 for fc € [1, s — 1] and Dgs = A. 
Hence the maximum of the sum ^k=ii^kk — 1) is Z\ — 1. 

The maximum of the product OLi i^k fe ~ 1) is reached when Dk fe = 2 for 
A: G [1, s — 1] and Dgs = so we have Y^k=ii^kk — 1) < — 1. 

The additional constraints on Dkk to be diagonal coefficients of the Smith 
normal form can only decrease an estimate. Next we trivially have n,. >oyk < 
nfe=i('^fefe ~ !)■ 

Moreover finding the Smith normal form U is a polynomial problem |d4j . 

Now we can prove the main result of this section. 

Theorem 7 Let A G b G P = P(A,b) be a simplex. If 

w{P{A,b)) > (5(A) — 1, then P{A,b) n Z" ^ 0. There is polynomial time 
algorithm to find some integer point in P. 

Proof Suppose that (5(A) > 1 (the case of d(A) = 1 is trivial). Let Ax < 6 be 
the subsystem of the system Ax < b, where A G and b G If such that 

|c?et(A)| = (5(A). Let C = {x G R" : Ax < b}. Then P c C. Let v G Q" be 
the vertex of P such that Av = b, so v = A~^b. 

Let B = —{6{A) — 1)A“^, then from elementary theory of polyhedrons it 
follows that C = v + cone{B). So columns of the matrix B defines edges of the 
cone C. Let S = v + conv{0 B). So S' C C. 

Also from elementary theory of polyhedrons it follows that n + 1 lines of 
the system Ax < b correspond to n + 1 facets of P by the following way: if F 
is a facet of P then F = P fl {a: G M" : Ak*x = bk} for some line (A^* bk) of 
the system Aa; < b. 

Let Fi, F 2 , Fn he the faces of P corresponding to the system Aa; < b 
and let P„+i be the last facet of P corresponding to the line of the system 
Ax < b that is not included to the system Aa; < b. 

First we need to prove that S fl Z" 0. Set C f If is induced by the 
following equivalent systems: 
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Ax < b 


{ Ax + y = b 
a; G Z" 
yGZ!^ 


(3) 


By the lemma U there is polynomial time algorithm to find a solution of this 
system y* G Z” with the property that X]fc=i Vk — '^(^) “ 1- Hence 3a:* G Z” 


that Ax* + y* = b and x* = A ^b — A ^y*. Finally a:* 


B 


1 


5(A) - 1 


y*, it 


is equivalent to the statement that a:* G S. 

To finish we need to prove that S C P. Suppose that S % P. Let , ..., 

be the vertices of S that are adjacent to the vertex v.li S % P then the facet 
Fn+i intersects some edge [a:, G S such that [u, fl Fn+i = {u}, 
u G vert{P) and u ^ Let Fj, where 1 < j < u, be an opposite facet 
to the vertex u and (Aj* bj) be the line of an initial system correspondent 
to Fj. Then |maa:{Aj*a: : x G P} — min{Aj^,x : x G P}| = \bj — Aj*u| < 

\bj — Aj I = 5(A) — 1. The last statement contradicts to an assumption of 
Theorem that w{P) > 5(A) — 1. 


Corollary 4 Let A G b G Z^^^, c P = P{A, b) be a simplex. 

Let the vertex v be an optimal solution of the linear problem max{c^x : x G 
P}, J = {j '■ Aj^,v = bj} and A = |det(Aj*)|. If w{P{A,b)) > A — 1, then 
P{A, 6) n Z" 7 ^ 0 and there is an algorithm with the eomputational eomplexity 
0(nA) that solves an integer problem max{c^x : a: G P fl Z"}. 


Proof We have already proved that S C P, where S = v + conv{{0 B)) 
and B = —{A — 1)A“^ (see the previous Theorem). It is easy to see that 
max{c^x : X G P n Z"} = max{c^x : x G S D Z"}. Also from Theorem [5] it 
follows that all vertices of the system @ are in the set S. But we already have 
a remarkable algorithm proposed by Gomory, Hu [Sl fTTlITB] for integer opti¬ 
mization in systems of the type o or m with a computational complexity 
0{nA). 

Corollary 5 Let A G 5 g Z^+^j c G Z" P = P(A, b) be a simplex. 

Ifw{P{A, b)) > Z\(A) —1, then P{A, 6)nZ"' ^ 0 and there is an algorithm with 
the computational complexity 0{nA) that solves an integer problem max{c^x : 
X G PnZ"}. 


4 Integer Programs with Almost Unimodular Matrices 

We will use the theorem [T] to prove the following result. 

Theorem 8 Let A G be almost unimodular, b G Z™, c G Z". Then, 

there is a polynomial-time algorithm to solve max{c^x : x G P{A, b) n Z"}. 
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Proof Let P = P{A, b). We can check an emptiness of P PI Z" in polynomial 
time by the theorem[TJ Then, if PflZ” ^ 0, we need to find a vertex v of P, that 
is an optimal solution of the relaxed problem max{c^x : Ax < b, x G Q”}. 
To this end, one can use the Khachiyan’s polynomial algorithm for linear 
programming |19j . If is integral, then it is an optimal solution of the integer 
problem. If not, we consider the shifted cone N{v) = {x € R” : Aj^^x < b}, 
where J = {k : A^^v = b}. 

By the theorem [TJ max{c^x : Ax < b, x G Z"} = max{c^x : x G 
N{v), X gZ'^} and an optimal integral point v* lies on some edge of N{v). 

Let Pfc be the set N{v) fl {x S R" : Xk = [I'fcl} and ^ be the set 
N{v) n {x G R" : Xk = bfcj}- 

Now we going to prove that there exists k, 1 < k < n such that v* is 
the optimal solution of the linear problem max{c^x : x G Bk} or problem 
max{c^x : x G Bk]. 

Let L be an edge of N{v). If |P fl Pfc| = 1 or |P fl B^ = 1 for some k and 
we know that LHBk = u or LCiBk = u then u G IP. Indeed, assume that L is 
induced by an integer subsystem of the system Ax < b denoted by Ax < b. We 
also can assume that A is {n — 1) x n matrix. Then the point u is induced by 
A \ /'b\ 

where tr is \vk\ or [ufc J. By assumptions 


the system ^ 

of a theorem the matrix 


is unimodular, hence u is integer. 


qA;—1 ^ Qn — k 

Hence we have only two possibilities for vert{Bk) (the same is true for 
vertjBk) ): 

1. vert{Bk) = {u}. This is the case when u G {a: G R" : Xk = so no 

new vertices were created. 

2. vert{Bk) C Z". Indeed if u ^ {a: G R" : Xk = then new vertices 

must be generated. Each vertex of Pfe is a result of an intersection of some 
edge of N{v) with the hyperplane {a; G R" : Xk = By the previous, all 

these intersections are integer. 

Now let L* be an edge of N{v) that contains an optimal point v*. 

There are few possible cases for some fixed k: 

1. L* C Bk. Then it is easy to see that an optimal solution of the linear 
problem max{c^x : x G Bk] is a point v. 

2. L* n Bk = 0. In this case the linear problem max{c^x : x G Bk] can 
be inconsistent or it has an integer or a rational solution. By the previous, if 
a solution is rational then it is the vertex v. 

3. |L*nPfc I = 1- This is the case of an intersection of a ray and a hyperplane. 
More precisely we can prove that L* Bk = {u*}. Let L* Ci Bk = {u} and 
u ^ V* . By the previous results the point u must be integral. Hence v* G [u, u] 
because L* is an edge of N{v) and an objective function only can be increased 
by the edge L*. But [u, u] can’t have any other integer points except u because 
Uk = r^fcl or Uk = [vfcj- So u = V*. 

4. The same items are true for L* Ci Bk. 
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Since L* contains v* G Z" then 3k, 1 < k < n such that \L* n Bk\ = 1 or 
\L* n = 1. Hence, by the previous, we have L* DB^ = {u*} or L* n ^ = 
{u*}. If this is true, then u ^ {x G M" : Xk = [ufc]} U {x G R" : Xk = [ufej} 
and all vertices oi Bk or B^ resp. are integer. 

In conclusion, we note that our final algorithm consists one running of 
a linear programming algorithm to find the vertex v as an optimal solution 
of the linear problem maxjc^x : Ax < b} and 2n times running of a linear 
programming algorithm to solve problems max{c^x : x G B^} and max{c^x : 
X G Bk} for each 1 < k < n. One of this solutions must be v* that can be 
recognized by an integrality and a maximality of the objective function. 

Since the linear programming is polynomial due to Khachiyan m so our 
algorithm is polynomial too. 

Let us to make some generalization of this proof. The matrix A is called k- 
almost unimodular for 1 < fc < n—1 if A{A) < 2 and any (r{A) — k) x {r(A) — k) 
sub-matrix has a determinant from the set {0, ±1}. 

Theorem 9 Let A G be /c-almost unimodular, b G Z™, c G Z". Then 

there is a polynomial-time algorithm for fixed k to solve max{c^x : x G 

P(H,&) nZ"}. 

Proof To proof this theorem we need to consider sets Bj = N{v) H {x G M" : 
xj = [uj]} and Pj = N{v) fl {x G K” : xj = [ujJ}, for all subsets J C I,n, 

I J| = k. There are (^) of these subsets. 

Again we have following possibilities for Bj: 

1. vert(Bj) = {u}. 

2. V* ^ vert{Bj) C Z". 

3. V* G vert{Bj) C Z". 

4. The same items are true for Pj. 

Since v* G Z” then 3J such that v* G vert{Bj) or v* G vert{^j). If this 
is true, then u ^ {x G M" : xj = [uj]} U {x G K” : xj = [ujJ}, hence the 
linear programs max{c^x : x G Bj} or max{c^x : x G Pj} have v* as an 
optimal solution. 

Our final algorithm consists one running of a linear programming algorithm 
to find the vertex v and 2 times running of a linear programming algorithm 
to solve problems max{c^x : x G Bj} and max{c^x : x G Bj} for each 
J C l,n, |J| = k. One of this solutions must be v* that can be recognized by 
an integrality and a maximality of the objective function. 


5 Examples of the Corner Polyhedrons with an Exponential 
Number of Edges 

Theorem [T] implies that there is a polynomial-time algorithm for solving an 
integer program max{c^x : x G P{A, 6)nZ"}, when A{A) < 2 and the number 
of edges in shifted cone N{v) is bounded by a polynomial, where v is optimal 
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solution of relaxed linear problem. But, there is example of an exponential 
number of edges in the shifted cone P{A, b), A{A) < 2. 

Consider the set D = conv{{0, 1}”} = {cc S K" : 0 < 2 :^ < 1, i € 1, n}. 

We use the homogenization and find the matrix Ad, such that: : 

< 0} = ■ X G D, xo > 0}. Finally, we create the matrix A'^^ by 

multiplying any column of Ad by 2. Ad is unimodular so A'jj is bimodular 
and we have the following theorem: 

Theorem 10 The cone {x G : A'j^x < 0} has an exponential number of 

edges. Any of its edges can be represented by a bimodular submatrix of Ad 
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